Accepted eor publication in The Astrophysical Journal Letters 
Preprint typeset using IATjiX style emulateapj v. 04/03/99 



RECENTLY DISCOVERED PULSARS AND UNIDENTIFIED EGRET SOURCES 
Diego F. Torres|, Yousaf M. Butt|, and Fernando Camilo| 

"Physics Department, Princeton University, NJ 08544 

*Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138 
^Columbia Astrophysics Laboratory, Columbia University, 550 West 120th Street, New York, NY 10027 
Accepted for publication in The Astrophysical Journal Letters 



o 
o 

(N 

a* 

C/3 



> 

oo 
(N 
(N 

O 



6 ■ 

a: 



ABSTRACT 

We present a correlative study between all unidentified EGRET sources at low Galactic latitudes and the newly 
discovered pulsars in the released portion of the Parkes multibeam radio survey. We note 14 positional coinci- 
dences: eight of these are "Vela-like" pulsars, with relatively small periods, small characteristic ages, and high 
spin-down luminosities. Three of these coincidences have been investigated by D' Amico et al. (2001) and Camilo 
et al. (2001). Among the others, we argue that PSR J1015— 5719 may plausibly generate part of the high energy 
radiation observed from 3EG J1014-5705. Three additional interesting cases are: 3EG J1410-6147 and either of 
PSRs J1412-6145 or J1413— 6141, if the pulsars are at the estimated distance of the coincident SNR G312.4-0.4; 
and 3EG J 1 639-47 02/PSR J1637-4642. The remaining positional coincidences between the EGRET sources and 
the newly discovered pulsars are almost certainly spurious. 

Subject headings: pulsars: general — pulsars: individual (PSR J1015— 5719, PSR J1412— 6145, 
PSR J1413-6141, PSR J1637-4642) — gamma rays: observations 



1. INTRODUCTION 

The Third EGRET Catalog lists 271 point sources (Hartman 
et al. 1999) of which about two thirds are not yet identified. 
Although most known young pulsars are not detected in the 
high-energy window (> 100 MeV), there are so far six such 
confirmed 7-ray pulsars (Thompson 1999). They are relatively 
young objects, with hard high-energy spectral indices, and ap- 
parently showing a trend for spectral hardening with increas- 
ing characteristic age (Fierro et al. 1993). Young pulsars were 
proposed as potential counterparts for 7-ray sources as early 
as 1983 (D' Amico 1983; Helfand 1994), and their relationship 
with unidentified EGRET sources remains a matter of current 
debate (see, e.g., Zhang, Zhang, & Cheng 2000). 

In order to unambiguously identify a pulsar as the origin of 
the 7-rays from an EGRET source, 7-ray pulsations must be de- 
tected at the pulsar period. However the 7-fluxes are generally 
very low, and nearly contemporaneous radio/X-ray and 7-ray 
observations are required to fold the few available photons with 
the correct ephemeris. Alas this is no longer possible owing to 
the demise of the Compton Gamma-Ray Observatory, and new 
candidate associations must for now be judged by comparison 
with the properties of the six known EGRET pulsars. 

Supernova remnants (SNRs) interacting with nearby molec- 
ular clouds have also been proposed as counterparts to some 
7-ray sources (e.g. Stumer, Dermer, & Mattox 1996; Esposito 
et al. 1996; Combi et al. 1998a, b, 2001), and in such cases 
both a pulsar and an associated interacting SNR may partially 
contribute to generate a given EGRET source flux. 

The Parkes multibeam pulsar survey is a large-scale sur- 
vey of a narrow strip of the inner Galactic plane (\b\ < 5°, 
260° < / < 50°; Manchester et al. 2001). It has much greater 
sensitivity than any previous survey to young and distant pul- 
sars along the Galactic plane, and it has resulted in the detection 
of many previously unknown young pulsars, potentially coun- 
terparts of EGRET sources. In this Letter we correlate the re- 



leased portion of the Parkes pulsar survey 1 , now consisting of 
368 pulsars, with all unidentified EGRET sources at low lati- 
tudes, in order to determine whether there are any new likely 
physical associations. 

2. CORRELATIONS 

Table 1 lists some important properties for the known 
EGRET 7-ray pulsars. Two of the parameters listed are highly 
uncertain: the distance d, and the observed efficiency for con- 
version of spin-down luminosity E into 7-ray luminosity L 7 : 
rj = L^/E = f4ird 2 F~ ( /E, where F 1 is the observed 7-ray flux 
and / is the 7-ray beaming fraction (0 < / < 1). Because this 
fraction is essentially unknown (see, e.g., Yadigaroglu & Ro- 
mani 1995), and in part because / ~ 1 results in -q ~ 100% for 
at least one pulsar at its nominal distance, it is common prac- 
tice, which we follow, to assume f = 1 /4ir. The efficiencies 
observed span the range 0.01% <r)<^ 3/19%: the two upper 
limits for 77 stem from the different measurements of the dis- 
tance to PSR B 1055-52 (see Ogelman & Finley 1993; Combi 
et al. 1997; Romero 1998). We shall then consider 0.01% £ 
77 <J 10% as "reasonable" values for evaluating the plausibil- 
ity of new associations, by analogy with presently known 7-ray 
pulsars, even though this inferred range relies on an assumed 
uniform value for /, and uncertain distances. 

All positional coincidences between unidentified 3EG 
sources and Parkes pulsars are given in Table 2. We list for each 
3EG source the variability index /: pulsars are steady 7-ray 
sources. For each pulsar we also present the "spin-down flux" 
E/d 2 , which is generally well-correlated with the detectability 
of 7-ray pulsars. A strict (linear) reading of this correlation 
would suggest L 7 oc E, while Thompson et al. (1999) find that 
L 7 oc E 1 / 2 . However, we summarily dismiss some potential as- 
sociations below owing to a low value ofE/d 2 only in extreme 
cases where the difference between E- and £' 1 ' 2 -scaUng is ir- 
relevant (and where even a hefty error in a more relevant and 



1 See http://www.atnf.csiro.au/research/pulsar/pmsur\ 
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uncertain parameter — distance — does not change matters 
significantly). In all these discarded cases, in addition, the ef- 
ficiencies required would be 77 3> 1000%, making the potential 
associations utterly unphysical. Finally, in evaluating possible 
associations, we compare the photon indices of known pulsars 
(Table 1) with those of the EGRET sources (see, e.g., Merck et 
al. 1996; Zhang & Cheng 1998; Cheng & Zhang 1998). 

We have also quantified the chance probability for obtain- 
ing these associations, adapting the numerical code described 
by Romero et al. (1999a, b): there is a 2% probability of hav- 
ing 8 chance coincidences between different unidentified 3EG 
sources and Parkes pulsars, as in Table 2. 2 

Camilo et al. (2001) proposed the possible physical associa- 
tion #1, with required efficiency of 0.5% (Table 2). This pulsar 
appears to be located just outside the SNR G284.3-1.8, which 
itself is interacting with an adjacent molecular cloud (Ruiz & 
May 1986). Another Parkes pulsar, PSR J1013-5934 (#2 in 
Table 2), is coincident with 3EG 1013-5915, but it is an old 
pulsar (r = 12Myr) with low E = 2.5 x 10 32 ergs s" 1 and cannot 
be a significant 7-ray contributor. D'Amico et al. (2001) have 
studied cases #6 and 14 (efficiencies: r] = 2% and 7%, respec- 
tively). Both are plausible candidates to generate the respec- 
tive EGRET source fluxes. PSR J 1837-0559, is also coincident 
with 3 EG J 1837-0606 (#13), but its E/d 2 is 80 times smaller 
than for PSR J 1837-0604, and it cannot contribute significantly 
to the 7-ray source. The pulsars in pairs #8 and 9 have far too 
low a spin-down luminosity at too large a distance (Table 2) 
to explain their coincident 7-ray sources. They are also both 
old, with r ~ 4Myr. Cases #10 and 1 1 are discussed elsewhere 
in connection with a proposal for a SNR shock origin of the 
bulk of the 7-rays resulting from 3 EG J1714-3857 (Butt et al. 
2001); neither of these pulsars is energetic enough to contribute 
significant amounts of high-energy flux. 

3. DISCUSSION 

We now discuss the remaining four EGRET sources position- 
ally superposed with five newly discovered pulsars. We provide 
observational data on the apparent associations in Tables 3 and 
4. Cases #3-5, 7 and 12 all contain Vela-like pulsars, with rel- 
atively short periods, low characteristic ages, and high spin- 
down luminosities (E ;> 10 35 ergs s" 1 ). 

The pulsar in case #3 would require an efficiency 77 = 5% 
at its nominal distance to explain the luminosity of the corre- 
sponding 3EG source, which has photon index 2.23 (see Ta- 
ble 4). This spectrum is softer than that of the Crab, although it 
is consistent with it within the uncertainties. It is also consistent 
with the index of 3EG 2227+6122, for which PSR J2229+61 14 
has been proposed as the likely source (Halpern et al. 2001). 
The 7-ray source in case #3 is not variable (Tompkins 1999; 
Torres et al. 2001c), as expected from direct pulsar or pul- 
sar wind nebula/SNR shock emission. PSR J 1015-57 19 has 
t = 39kyr and £ = 8.2 x 10 35 ergs s" 1 (Table 3). While no cat- 
aloged SNR is superposed with the 3EG source (Torres et al. 
2001b), this absence does not mean that one does not exist, and 
further sensitive searches may prove fruitful. Thus the connec- 
tion between 3EG J1014-5705 and PSR J1015-5719 appears 
plausible and is worth additional study. 

The pulsars in pairs #4 and 5 require unreasonably high ef- 
ficiencies at their nominal distances to explain the 7-ray flux 
from the corresponding EGRET source (77 ^ 100%, Table 4). 

2 To arrive at this result we have considered the 41 3EG sources in the range of 
average expected number of chance coincidences. There is therefore a 98% a prio 



However, both pulsars are located in the direction of the Cen- 
taurus arm, and it is known that in such directions the elec- 
tron density/distance model of Taylor & Cordes (1993) can 
be unreliable, sometimes overestimating the distances by fac- 
tors of up to <~ 4 (see discussion in Camilo et al. 2001). 
Both pulsars are located (at least in projection) well within the 
boundaries of the incomplete shell SNR G3 12.4-0.4 (Caswell 
& Barnes 1985), to which Yadigaroglu & Romani (1997) es- 
timate aS-D distance of 1.9 kpc. At this distance the re- 
quired efficiencies for the pulsars in cases #4 and 5 would 
be 12% and 3%, respectively, which would make them con- 
siderably more plausible sources of the observed high energy 
emission. Furthermore, 3EG J1410-6147 has a photon index 
comparable to that of the Crab, and is not variable (Tables 1, 
2 and 4). Pairs #4 and 5 therefore appear intriguing. How- 
ever it should be noted that E-D distances are notoriously 
unreliable (e.g., for this very SNR, Caswell & Barnes 1985 
and Case & Bhattacharya 1999 infer values in substantial dis- 
agreement both with each other and with that determined by 
Yadigaroglu & Romani 1997). Ideally, further observations of 
SNR G3 12.4-0.4 may indicate whether it shows signs of in- 
teraction with PSRs J1412-6145 or J1413-6141, and possibly 
constrain their distances. Depending on the actual distances, 
the 7-ray emission from 3EG J1410-6147 may conceivably 
arise from a combination of PSRs J1412-6145, J1413-6141, 
and/or SNR G3 12.4-0.4. 

The efficiency required to explain the EGRET flux in case #7 
is 77 = 12%, which seems possible. However, the spectral index 
of 2.50 is larger than those of known 7-ray pulsars. One of the 
SNRs coincident with the EGRET source, G337.8-0.1, harbors 
a maser (Koralesky et al. 1998), which is indicative of interac- 
tion between the SNR shock and the ambient medium. Thus, 
were a sufficiently massive molecular cloud located nearby, 
it could help produce the high energy radiation as a result of 
hadronic interaction (Aharonian, Drury, & Volk 1994; Aharo- 
nian & Atoyan 1996). Part of the EGRET flux could plausi- 
bly come from PSR J 1637-4642 and part from pion 7-decay 
via SNR G337.8-0.1's interaction with the putative cloud. In 
this case, the photon index would reflect a weighted average 
value. Both possible mechanisms for the high energy emis- 
sion would produce a non-variable source, as is the case for 
3EGJ1639-4702. 

Lastly we consider case #12, for which the required effi- 
ciency is high at the nominal pulsar distance and upper limit 
flux value, i] = 55%. The pulsar is located, in projection, just 
outside the plerionic SNR G27.8+0.6, for which the distance 
is ~ 2 kpc (Reich et al. 1984). Although the estimated ages 
are comparable (t = 52kyr for the pulsar and ~ 45kyr for the 
SNR), it seems unlikely that both objects are physically associ- 
ated, given the offset between the centrally peaked SNR com- 
ponent and the pulsar (see Reich et al. 1984). Whatever the 
possible relation between pulsar and SNR, the 3EG source in 
case #12 is variable (Table 2), arguing against a pulsar origin. 

Examination of X-ray archives via HEASARC has revealed 
no compelling counterpart sources to any of the pulsars in Table 
3. (A possible Einstein X-ray source at the edge of an IPC field 
near the location of PSR J1015-5719 [#3] is likely an artifact 
resulting from the large effective area correction used.) How- 
ever, the exposures may not have been deep enough to reveal 
such X-ray emission, usually expected at the level L x ~ 10~ 3 E 
(e.g., Becker & Trumper 1997). Future observations of these 

he Parkes survey, and simulated them in synthetic populations in order to obtain an 
i probability that not all 8 associations are spurious. 
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pulsar fields with XMM-Newton and Chandra could thus be in- 
structive. 

In conclusion, we find that two recently discovered Parkes 
pulsars (PSRs J 1015-57 19 and J 1637-4642) could plau- 
sibly generate at least part of the 7-ray flux observed 
from two unidentified EGRET sources (3EG J1014-5705 
and 3EG J 1639-4702, respectively), and either of 
PSRs J1412-6145 or J1413-6141 could be associated with 
3EG J1410-6147 if they are located closer than their DM dis- 
tances by a factor of ~ 4. Cases #3-5 and 7 represent promising 



targets for the forthcoming AGILE and GLAST missions. 
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Table 1 

7-ray pulsars detected by egret 



Pulsar/3EG Jsource 


P 

(ms) 


T 

(kyr) 


E 

(ergss -1 ) 


d 
(kpc) 


F^ EG [xl0- s ] 
(ph cm" 2 s" 1 ) 




V 

(100MeV- 
-lOGeV) 


Crab/0534-2200 


33 


1.2 


5.0 x 10 38 


2.0 


226.2±4.7 


2.19±0.02 


0.01% 


Vela/0834-4511 


89 


12.5 


6.3 x 10 36 


0.25 


834.3±11.2 


1.69±0.01 


0.08% 


B1951+32/ ••• 


39 


100.0 


3.7 x 10 36 


2.4 






0.3% 


B1706-44/1710-4439 


102 


15.8 


3.1 x 10 36 


1.8 


111.2±6.2 


1.86±0.04 


1% 


Geminga/0633+1751 


237 


316.2 


3.1 x 10 34 


0.16 


352.9±5.7 


1.66±0.01 


3% 


B1055-52/1058-5234 


197 


501.1 


3.1 x 10 34 


0.5/1.5 


33.3±3.8 


1.94±0.10 


2/19% 



Note. — Pulsar parameters and distances (and rj in the case of PSR B 195 1+32 — this pulsar is not a 3EG 
source) are taken from Kaspi et al. (2000), excepting PSR B 1055-52, for which we also consider a smaller 
value of distance (Ogelman & Finley 1993; Combi et al. 1997), and Vela (Caraveo et al. 2001 and references 
therein), r = P/2P, and E = 4ir z IP/P\ with / = 10 45 gem . The "P1234" 7-ray fluxes and spectral indices 
are from the 3EG catalog (Hartman et al. 1999), from which we have computed the r\ values. 
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Table 2 

Positional coincidences between unidentified 3EG sources and Parkes pulsars 



JCU J 


iNote 


i 
i 


ncn I 
roK J 


if 
ft 


vMNKS 


Ker. 


Z7 / j2 

ti / a 


1013-5915 


C, em 


1.6 


1016-5857 


1 


G284.3-1.8 


Camilo et al. (2001) 


2.9 x 10 35 








1013-5934 


2 






2.0 x 10 30 


1014-5705 


C, em 


1.4 


1015-5719 


3 






3.4 x 10 34 


1410-6147 


C 


1.2 


1412-6145 


4 


G3 12.4-0.4 




1.4 x 10 33 








1413-6141 


5 






1.6 x 10 33 


1420-6038 


C 


2.1 


1420-6048 


6 




D'Amicoetal. (2001) 


1.7 x 10 35 


1639-4702 


C, em 


1.9 


1637-4642 


7 


G337.8-0.1; G338. 1+0.4; G338.3+0.0 




1.9 x 10 34 








1640-4648 


8 






1.5 x 10 32 








1637-4721 


9 






3.1 x 10 30 


1714-3857 


C, em 


2.1 


1713-3844 


10 


G348.5+0.0; G348.5+0.1; G347.3-0.5 


Butt et al. (2001) 


4.0 x 10 31 








1715-3903 


11 




Butt et al. (2001) 


3.0 x 10 33 


1837-0423 


C 


5.4 


1838-0453 


12 


G27.8+0.6 




1.2 x 10 33 


1837-0606 


C, em 


2.4 


1837-0559 


13 






6.5 x 10 32 








1837-0604 


14 




D'Amico et al. (2001) 


5.2 x 10 34 



Note. — A number (#) is assigned to each pair for ease of reference. SNRs contained in Green's (2000) catalog coinciding with 
the EGRET sources are noted (see Torres et al. 2001b). Only in cases #4, 5 and 12 do the pulsars also coincide with the SNRs 
listed. "C" and "em" refer to the 7-ray sources: source confusion exists and sources are possibly extended or multiple, respectively 
(Hartman et al. 1999). / is the variability index as in Torres et al. (2001a, c), where / > 5 (< 2) represents a source whose flux 
presents variability levels at least 8 a (less than 2a) above those displayed by confirmed pulsars. Pulsar "spin-down flux" E/d 2 is 
in units of ergs s" 1 kpc" 2 . 



Table 3 

Observational parameters for the pulsars in the candidate associations 



Case 


A9 [deg] 


[deg] 


(l,b) 


d [kpc] 


r [kyr] 


F[ms] 


P [10" 15 ] 


E [ergss '] 


#3 


0.30 


0.67 


283.09,-0.58 


4.9 


38.7 


140 


57.4 


8.2 x 10 35 


#4 


0.15 


0.36 


312.32,-0.37 


9.3 


50.6 


315 


98.7 


1.2 x 10 35 


#5 


0.28 


0.36 


312.46,-0.34 


11.0 


13.5 


286 


333.4 


5.7 x 10 35 


#7 


0.46 


0.56 


337.79,+0.31 


5.8 


41.2 


154 


59.2 


6.4 x 10 35 


#12 


0.50 


0.52 


27.07,+0.71 


8.2 


52.2 


381 


115.7 


8.3 x 10 34 



Note. — A8 is the angular distance between the center of the 3EG source and the position of the re- 
spective coincident pulsar. 9 is the effective 95% confidence level radius of the 3EG source error box 
(Hartman et al. 1999). Pulsar parameters are taken from the Parkes survey database (see footnote 1), and 
their distances (d) are estimated from the observed dispersion measure (Taylor & Cordes 1993). 
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Table 4 

Observed and computed 7-ray emission parameters for new candidate associations 



Case 


F 3EG [xlO _8 ph cm" 2 s" 1 ] 


p3EG [ er g Scm 2 S '] 


7 3£G 


^isotropic [^*gS S ] 


Lbeamed [ergs S '] 


7] 


#3 


34.0±6.5 


1.90 xlO" 10 


2.23±0.20 


5.2 xlO 35 


4.1 xlO 34 


5% 


#4 


64.2±8.8 


4.09 xlO" 10 


2.12±0.14 


4.1 xlO 36 


3.2 xlO 35 


270/1 2% a 


#5 


64.2±8.8 


4.09 xlO" 10 


2.12±0.14 


5.7 xlO 36 


4.5 xlO 35 


80/3 % a 


#7 


53.2±8.7 


2.30 xlO" 10 


2.50±0.18 


9.3 xlO 35 


7.4 xlO 34 


12% 


#12 


<19.1 


< 0.70 x 10" 10 


2.71 ±0.44 


<5.8 x 10 35 


<4.6x 10 34 


<55% 



Note. — The EGRET fluxes correspond to the "P1234" values (Hartman et al. 1999). Luminosities, and efficiency values for 
all pulsars, are given for the distance estimated from the dispersion measure (Table 3). 

a Second values of efficiency correspond to a distance of 1.9kpc, one of those estimated for the positionally coincident 
SNR G3 12.4-0.4 (see text). 



